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CHAPTER 1 - INTRODUCTION 


1.1 An Overview of the Thesis 

The problem of determining a maximum steady state flow from one 
point to another in a network with edge capacity constraints, has come to 
be known as the maximum network flow problem (“max-flow" in short). L.R. 
Ford and D.R. Fulkerson [9] were the first to study max-flow as a 
computational problem. They developed the first max-flow algorithm in the 
mid 1958’s and laid the groundwork for much of the research that was to 
follow. Since that time max-flow has been widely studied and has developed. 
great practical application, especially in the analysis of transportation 
and communication networks. This thesis will be concerned with the 
computational complexity of the maximum network flow problem, investigating 
both upper and lower bounds on the problem. 

The remainder of this chapter will be devoted to a development of 
the foundation necessary for any coherent study of the maximum network flow 
problem. We will begin by presenting the basic definitions and concepts 
that have become standard in the max-flow literature. We shall then 
introduce the fundamenta! theorems upon which the study of max-flow has 
been built. These theorems, due originally to Ford and Fulkerson [9], are 
known as the Augmenting Path Theorem, the Integral Flow Theorem and the 
Max-Flow Min-Cut Theorem. 

The second chapter in this thesis will deal with upper bounds on 
the complexity of the maximum network flon problem. The long and 


intriguing history of the search for such bounds is summarized below in 


Table 1.1. 
Upper Bounds on Max-F low 
1) Ford and Fulkerson (1956) Unbounded 
2) Edmonds and Karp (1969) O(V-E2) 
3) Dinic (1978) ov) 
4) Karzanov (1973) _ o(v3) 
5) Cherkasky (1976) ov2._l/2) 
6) Galil (1978) o(v9/3.¢2/3) 
7) Malhotra, Kumar and Maheshwari (1978) o(Vv3) 
8) Galil and Naamad (1978) O(V-E-1og2V) 


Table 1.1 


Each of these bounds has been demonstrated by the construction of a max- 
flow algorithm which has the specified worst case running time. Chapter 2 
will consist of a survey of each of these algorithms. We note here that 
the algorithm developed by Malhotra, Kumar and Maheshwari [16], which we 
will call the MKM algorithm, does not result in an asymptotic improvement 
over the previous three algorithms. However, this algorithm is very simple 
and is probably the best algorithm to use on dense netunorks (ExV*), 

The determination of lower bounds on the computational complexity 
of the maximum netuork flow peentew has thus far received little attention 
in the Livenatucs: In Chapter 3, however, we shall investigate one 
particular approach to establishing a non-linear lower bound on the 
complexity of max-flow. The technique we shall deal with is the polyhedral 
technique developed in 1977 by A.C. Yao, O.M. Avis and R.L. Rivest [21]. 


We will show that a straightforward application of this technique fails to 


produce a non-linear {fower bound on max-flow. In the process, however, we 
shall answer several questions concerning the facial structure of a class 


of polyhedra very closely related to the max-flow problem. 


1.2 Basic Definitions and Concepts 
This section will be devoted to a presentation of the basic 


definitions and concepts that have become standard in the study of max- 


flow. We will begin by defining a network and a legal flow function on a 
network. These concepts will then be used to develop a, formal definition 
of the maximum network flow problem. Finally we will introduce the basic — 


notions of a cut and a flow augmenting path. 


Definition 1.1: 
A netunork = (G,s,t,c) is a 4-tuple with the properties: 


1) G = (V,E) is a finite directed graph composed of a 
set of vertices, V, and a set of edges, E. 


2) Two distinct vertices seV and teV are specified as 
the source and sink respectively. 


3) Each edge eeE is assigned a non-negative real number 
cle), called the capacity of edge e. 


For any vertex veV, In(v) denotes the set of all edges incoming to 
v and Out(v) denotes the set of all edges outgoing from v. An edge e, 
directed from a vertex u to a vertex v in a netnork M, will often be 


denoted by the ordered pair e=(u,v). 


Example 1.1: (Network Ng) 


V = {s,a,b,c,d, t} 

E = }(s,a),(s,b), (a,c), (a,d), 
(b,d), (c, t), (a, t)} 

cle) = 18 [for all ecE s.t. ex(a,d)] 

cla,d) = 1 

In(t) = {(c, t), (d, t)} 

Out(s) = {(s,a), (s,b)} 


Definition 1.2: 
A legal flow function f on a network 7 associates with each edge 
ee€E a real number fle) satisfying the conditions: 
C1) Osfle)scle), for each edge eeE 


C2) >» fle) - > fle) = 8, for each vertex veV-{s, t}. 
eelIn(v) eeQut (v) 


The value of a legal flow function f is defined to be: 


(1.1) vif) = D fle) - > fled. 
eeIn(t) eeQut (t) 


Informally we say f is a steady state flow from s to t and fle) is 
the steady state flow through edge e. Thus, condition Cl telis us that the 
steady state flow through any edge must be non-negative and not exceed the 
capacity of that edge. Further, condition C2 states that any steady state 
flow from s to t must have the property that: flow is conserved at every 


vertex other than the source and the sink. 


Definition 1.3: 
An edge e, incident upon vertices u and v in a network 7, is said 


to be useful from vertex u_to vertex v with respect to a legal flow 


function f defined on Ml, if either: 


1) e=fu,v) and fle)<cle) 
or 
2) e=(v,u) and fle)>g. 


Example 1.2: 


The following tegal flow function fg is defined on network Ng from 
Example 1.1: . 


18,4 
.fg(s,a) = 5 fg (s,b) = 8 
fgla,c) a4 fg la,d) = 1 
fg(b,d) = 8 fg{c, t) = 4 
fg (d, t) a j 


18,8 


v (fg) = fg{c, t) + fg(d, t) - 8 = 44] #5 
Note - (s,b) is useful from s to b but not useful from b to s 


(a,d) is useful from d to a but not useful from a to d 
(a,c) is useful from a to c and from c to a 


Let us now consider a legal flow function, defined on a network 7, 
whose value is maximum over the set of values of all legal flow functions 
defined on It. Such a flow function is said to be maximum with respect to 
Mm. Notice that the existence of a maximum flow function on a network is 
not open to question since a network is composed of only finite capacity 
edges. Further, notice that there may exist more than one distinct maximum 


flow function on a network. 


Definition 1.4:  (Max-Flow) 


The maximum network flow problem is defined as the problem of 
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computing a maximum flow ‘function on a network given as input. 


We now turn our attention to the basic concepts underlying the 
theory on which the study of the maximum network flow problem has been 
built. In particular, we present the notions of a cut and a flow 


augmenting path. 


Definition 1.5: 
A cut in a network It is a set of vertices X with the properties: 
1) XceV. 
2) sex, 
3) teX where K = V-X. 
The set of all edges eeE which are directed from a vertex in X to a vertex 
in X is denoted by (X;X). The capacity of a cut X is defined to be: 


(1.2) CX) = DS cle). 
ee (X;X) 


Example 1.3: 
The following cut Xg is defined on network Ng from Example 1.1: 
4 
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Xg = {s, a, d} 


Xg = {b,c, t} 


C{Xg) = c(s,b) + cla,c) + cld,t) = 18+18+19 = 30 


ii 


Definition 1.6: 


A path p, from vertex v, to vertex v. ina network M, is a 


n 
sequence of distinct vertices and edges p = Vj@)V5e9 «++ Vp__yOn_1 Yh, (n22) 
such that v eV (for each ixl,...,n), e,€E (for each ixl,...,n-1) and 
either e; = (VisViyg)s in which case e; is said to be a forward edge inp, 


or e; = (viygeVj)> in which case e; is said to be a reverse edge in p, for 
each i=l,...,n-l. Any path composed of only forward edges is called a 


chain. Note that a path is defined to be acyclic. 


Definition 1.7: 

A flow augmenting path, with respect to a legal flow function f on 
a network 7I?, is a path p” with the property that each edge e, ep” is useful 
from vertex v,ep* to vertex Vig} EPs A flow augmenting path from the 
source to the sink is called an s-t flow augmenting path. We shall see in 
the next section that an s-t flow augmenting path is a path along which f 


can be augmented (i.e. the value of f can be increased). 


Example 1.4: 
The following s-t flow augmenting path Pg is defined with respect 


to the legal flow function fg from Example 1.2: 
Pg = s(s,b)b(b,d)dla,d)ala,c)clc, t)t 


Note - (s,b), (b,d), (a,c) and (c,t) are all forward edges along Pg 
(a,d) is a reverse edge along Pg: 


12 


1.3 Fundamental Theorems 

The study of the maximum network flow problem has been based on 
three fundamental theorems first proven by Ford and Fulkerson [9] in the 
mid 1958's. These theorems are known as the Augmenting Path Theorem, the 
Integral Flow Theorem and the Max-Flow Min-Cut Theorem. This section will 
be used to develop formal proofs of each of these theorems. We will begin 
by proving the following principle Jemma which we alluded to at the end of 


the last section. 


Lemma 1.1: 
If there exists an s-t flow augmenting path p’, with respect to a 
legal flow function f on a network M, then there exists a legal flow 


function f% on Ml such that v(f“)ov(f). 


Proof: 
Let PY EV, C4 V9E9+ + Vn On Vy {n22) be an s-t flow augmenting path 
With respect to a legal flow function f on a network It. Define the sets 
Eg. Ey and E> as follows: 
Eg = fecE |e is not an edge in p%} 
E, = feeE |e is a forward edge in p*} 
E> = feeE |e is a reverse edge in p*}. 
Note that E,UE> is the set of all edges in p’ and that EglE, UES is the set 
of all edges in 7 (i.e. EguE,UE> = E). 
Since p’ is an s-t flow augmenting path, it follows from Definition 
1.7 that V1] =3, Vnzt and e; is useful from v. to Vial? for each i=#l1,...,n-l. 


More explicitly, applying Definition 1.3 yields: 
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(1.3) fle)<cle) = cle)-fle)>8, for each edge eeE, 
and 
(1.4) f(e)>@8, for each edge eeE5. 


Now construct the function f” as follous: 


(1.5) Cle)ec(e)-fle), for each eeE, 
(1.6) clelef(e), for each eee, 
& © min cle) 
e€E, UE, 
(1.7) f’(elefle), for each ecEg 
(1.8) f’(e)ef(el+8, for each eek, 
(1.9) f’(e)efle)-8, for each e€E>. 


It must now be shown that f” is a legal flow function on WM and that 
v(f“)>v(f). (Notice that f” is defined on all edges eeE.) 

To show that f% is a legal flow function on M, we must simply 
prove that f% satisfies both condition C1 and condition C2 of Definition 
1.2. Let us first consider condition Cl. Since f is a legal flow function 
on [t, it follous from Definition 1.2 that: 

(1.18) O<sfle)<cle), for ali eee. 
Thus, applying (1.7) yields: 
(1.11) 8st’ le)scle), for all ecEp. 
We now notice, from (1.3) and (1.4), that cle)>8 for all ecE,UE> and thus 
6>8. Combining this with (1.8) and (1.9), we have: 
(1.12) f’(e)>fle), for all eeE; 
and 


(1.13) f’le)<fle), for all eeE>. 
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We can now apply (1.18) and (1.12) to (1.8) to obtain: 
(1.14) @<f’le)<cle), for all eeE;. 
Similarly, we can apply (1.18) and (1.13) to (1.9) to obtain: 
(1.15) Osf’(el<cle), for all ecE>. 
Therefore, combining (1.11), (1.14) and (1.15), we have: 
Bs’ (e)scle), for all eeE 

and thus f” satisfies condition Cl of Definition 1.2. 

We must now prove that f” satisfies condition C2 of Definition 1.2. 
This is accomplished by considering any vertex v’eV-{s,t}. Since f is a 
legal flow function on 7, it follows from Definition 1.2 that: 


(1.16) > fle) - DY fled = 8. 
a eelIn(v’) eeQut (v”) 


Further, if v“ is not a vertex along the path p’, then all of the edges 
incident upon v” must be contained in the set Eg. Thus, combining (1.7) and 
{1.16) yields: 
> tte) - Rete) = Hetey - DY fled - 8. 
eelIn(v’) eeQut(v*) © eeIn(v’) ecOut (v") 


(for all v’eV-{s,t} and v“¢p%) 


If, however, v’ is a vertex along the path p*%, then let i be the index of 
v’ along p’ (i.e. v’ = v;ep’). Since p” is acyclic and v’¢{s,t}, there 
must be exactly two distinct edges incident upon v” which are contained in 
p’, namely €;_, and e;. All other edges incident upon v” are therefore 
contained in the set Eg: We must now consider each of the following four 


possible cases: 


, 


1S 


1) e;_,€E) and e;€E, 
2) e, 1} and e, €E> 
3) e;_,€E> and e; EF; 
4) e;_;€E> and e;€E>. 
For case 1, we combine (1.7), (1.8) and (1.16) to obtain: 


> fle) - DY tte -(& rier) + 6] - (> ted) + 8] 


eeIn(v’) eeDut (v7) eeIn(v% eeOut (v%) 
(Ste) - (Sete) + 8 - 
eeIn(v’) eeOut (v7) 


- Detter - DS tte) = 8, 


eeIn(v’)  eeDut(v’) 
The same result is proven for cases 2, 3 and 4 similarly. We therefore 
have that: 


> tite) - DS tle) = 8, for ali veV-{s, t} 
eel n(v”) eeDut (v7) 


and thus f” satisfies condition C2 of Definition 1.2. 

To show that vi(f’)ov(f), we first notice that there is exactly one 
edge incident upon t which is contained in the path p’, namely en-1° All 
other edges incident upon t are therefore contained in the set Eg: We must 
now consider each of the following two possible cases: 

1) @ 4.4 €E} 
2) e,_;6E2- 
For case 1, we can simply combine (1.1) of Definition 1.2, (1.7) and (1.8) 


to obtain: 


v(t) = D> tle) - DY t’ led 


eel n(t) eeQut (t) 


(> te) + 5] - D> tle) 


eelIn(t) eeQut (t) 


16 
~ Dette) - D tte) + 8 
eeIn(t) eeOut (t) 
= vif) +8 > vf). (since 8>8) 
The same result is proven for case 2 similarly. We therefore have that 


v(t“) > v(t). 


The proof. of Lemma 1.1 reveals a fairly simple procedure for 
augmenting an existing legal flow function on a network, given a 
corresponding s-t flow augmenting path. The procedure consists of 
determining the "excess capacity" along each edge in the flow augmenting 
path, as defined by (1.5) and (1.6), and then increasing the flow along all 
forward edges and decreasing the flow along all reverse edges in the path. 
The amount by which flow is increased or decreased along each edge is 
simply the minimum excess capacity over all edges in the flow augmenting 
path. 

We now present a lemma which will be useful in demonstrating a 
relationship between the capacity of a cut and the value of a legal flow 


function on a network. 


Lemma 1.2: 


Given any cut X and any legal flow function f on a network 7, 


vif) = DD tle) - DY fled. 


e€ (X;X) ee (X;X) 


Proof: 


Let X be any cut and f be any legal flow function on a network 7. 
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Then by condition C2 of Definition 1.2 and the fact that s ¢ X, we have: 


(1.17) > (S tte - Herter) = 2. 


veX-{t} eeln(v) eeOut (v) 


Combining (1.1) of Definition 1.2 and (1.17) yields: 


1.189 vith = Hele - Deter + DY (Ket - > tte). 


eelIn(t) eeDut (t) vex-{t} eeIn{v) eeDut (v) 


Simplifying (1.18) we obtain: 


vif = DY (Seer - J eter) 


vex eelIn(v) eeQut (v) 
= >> > f (e) - > > fle). 
veX eeln(v) vex eeQut(v) 


Thus, we have that the value of f is equal to the sum of the flow along all 
edges directed into a vertex in X minus the sum of the flow along all edges 
directed out of a vertex in K. If we now consider separately those edges 
directed from a vertex in X to a vertex in X and those edges directed from 


a vertex in X to a vertex in X, we obtain: 


vif) = (Setter + Het) - (Dee + ¥eter) 


ee (X;X) ee (X3X) ee (X;X) ee (K;X) 


= Dtted - Dd fled. 


ee (X3X) ee (X;X) 


Corollary 1.1: 
If the value of a legal flow function f” is equal to the capacity 
of some cut X”% on a network 7, then f” is maximum on I and X” has minimum 


capacity over all cuts on 7. 


Proof: 


Let X be any cut and f be any legal flow function on the network 
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r. Then by Lemma 1.2 we have that: 


vif) = > fled - DD fled. 


e€ (X3X) e€ (X;X) 
Applying condition Cl of Definition 1.2 we obtain: 


v(f) os > cfe) 
ee (X;X) 


and thus by (1.2) of Definition 1.5: 

(1.19) v(f) < CM). 

Therefore, if (1.19) holds by equality for some legal flow function f” and 
some cut X”% on the network M, then f”% must be maximum on Tt and X” must 


have minimum capacity over all cuts on 7, 


We are now ready to present the three fundamental theorems upon 


which the study of max-flow has been built. 


Theorem 1.1: (Augmenting Path Theorem) 
A legal flow function f on a network I? is maximum if and only if 


there exists no s-t flow augmenting path with respect to f on 7t. 


Clearly, if there exists an s-t flow augmenting path with respect 
to a legal flow function f on a network 7M, then by Lemma 1.1 f is not 
maximum on 7. Assume now that there is no s-t flow augmenting path with 


respect to f on It and define the set S as follows: 


S = {veV | J an s-v flow augmenting path with respect to f on 7} u {s}. 
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Since there is no s-t flow augmenting path with respect to f on 7%, it can 
easily be seen from Definition 1.5 that S forms a cut in 7. Further, from 
Definition 1.3 and Definition 1.7 we have that fle)=cle) for each edge 
e€(S;5) and fle) = 8 for each edge e€(5;S). We can now apply Lemma 1.2 to 


obtain: 


vit) = D fled - >} tled 


ee (5;5) ec (53S) 


= >: cle) 


e€(S;5) 
= C(S). 


Thus, by Corollary 1.2 we have that f is maximum on 7. 


Theorem 1.2: (Integral Flow Theorem) 
There exists an integral valued maximum flow function on any 


network defined by an integral valued capacity function. 


Let M be any network defined by an integral valued capacity 
function and Jet fg be the zero flow function on Il, defined by fg (e) =8 for 
each edge eeE (Notice that such a flow function will be a legal flow 
function on any network). We can now compute a maximum flow function on 7 
as follows: | 

while there exists an s-t flow augmenting path 
with respect to fy on Tt 
do augment fy as outlined in 


the proof of Lemma 1.1. 


28 


An examination of (1.5) through (1.9) reveals that the legal flow function 
generated at each iteration of this procedure will be integral valued. 
Thus by Lemma 1.1, the value of each successive legal flow function 
generated must be at least one integral unit greater than the value of the 
previous legal flow function. Combining this fact with (1.19) and Theorem 
1.1, we now have that our procedure must halt within a finite number of 


steps, yielding an integral valued maximum flow function on 7M. 


Theorem 1.3: (Max-Flouw Min-Cut Theorem) 
The value of any maximum flow function on a network 2 is equal to 


the minimum cut capacity over all cuts on 7. 


Let f be any maximum flow function defined on a network 7. 
Applying Theorem 1.1 we have that there is no s-t flow augmenting path with 
respect to f on 7. It now follows immediately from the proof of Theorem 
1.1 that there exists a cut S on Mf such that: 

v(f) = C(S). 
Further, by Corollary 1.2 we have that the cut S must have minimum capacity 
over all cuts on 7. Thus, the value of any maximum flow function on 7 is 


equal to the minimum cut capacity over all cuts on 7. 


We shall see in Chapter 2 how the previous three theorems form the 


basis for all the max-flon algorithms that have thus far been developed. 
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In fact, we will see that each algorithm is actually a variation of the 


procedure given in the proof of Theorem 1.2. 
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CHAPTER 2 - UPPER BOUNDS ON MAX-FLOW 


2.1 Introduction 

The formal definition of the maximum network flow problem, as 
presented in Section 1.2, can be explicitly stated as follows. Given any 
network I as input, compute values of the variables Xe [for each ee€E) so 


as to maximize the objective function 


y % - D Xe 


eelIn(t) eeQut (t) 
subject to the constraints 
X, 2 8, for each eeE 


X, < cle), for each eeE 


e 
>, Xe - > X, = 8, for each veV-{s, t}. 
eeln(v) eeDut (v) 


Thus the maximum network flow problem can be viewed as an optimization 
problem in which a linear function must be maximized subject to a system of 
linear equations and linear inequalities. G.B. Dantzig [4] developed an 
algorithm in the early 1958's, known as the simplex method, which could be 
used to solve such linear programming problems. Although it would not be 
incorrect to consider the simplex method to be the first max-flow 
algorithm, it is usually not treated as such. The simplex method is a very 
general algorithm which has an unbounded worst case running time. We have 
included it here only for the sake of completeness. 

L.R. Ford and D.R. Fulkerson (9) were the first to produce 
significant research results concerning, specifically, the maximum network 


flow problem. In the mid 1958's they proved the Augmenting Path Theorem, 
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the Integral Flow Theorem, and the Max-Flow Min-Cut Theorem. These 
fundamental results led directly to their development of the labeling 
algorithm for solving the maximum network flow problem. The labeling 
algorithm is a straightforward algorithm which simply augments an existing 
legal flow function along some s-t flow augmenting path in a network. This 
process is then repeated until there no longer exist any s-t flow 
augmenting paths in the network. The labeling algorithm, although it also 
has an unbounded worst case running time, remained in successful use for 
almost 15 years. 

In 1969, J. Edmonds and R.M. Karp [7] developed a variation of the 
labeling algorithm which utilized a Breadth First Search in picking out the 
s-t flow augmenting paths in ghdeeot increasing length. This resulted in 
a much more efficient algorithm with a bounded OV] -JEI2) worst case 
running time. Independently and a short time later, E.A. Dinic [5] 
developed an improved version of Edmonds and Karp’s algorithm. Dinic also 
utilized the technique of Breadth First Search but he developed an 
algorithm with time complexity OCIVIZ4E}). 

A.V. Karzanov [14] modified Dinic’s algorithm in 1973 to obtain an 
ocvy3) max-flow algorithm. Karzanov’s algorithm was unique in that it 
simultaneously augmented an existing legal flow function along several s-t 
flow augmenting paths. In 1976, B.V. Cherkasky [18] showed how to combine 
Dinic’s algorithm with Karzanov’s algorithm to produce a new and very 
complex Oty} 2-yeyt/2) max-flow algorithm. Two years later, Zvi Galil 18) 
improved Cherkasky'’s algorithm to ocvyo/3.jey2/3) by developing a 


technique for retaining useful information about the structure of the 
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network. 

In 1978 V.M. Malhotra, M. Pramodh Kumar and S.N. Maheshwari [16] 
discovered a very simple ovis) max-flow algorithm which we shall call the 
MKM algorithm. Their algorithm, similar to Karzanov’s algorithm in that it 
simultaneously augments along several s-t flow augmenting paths, requires 
very little poeanead: Although the MKM algorithm does not result in an 
asymptotic improvement over the previous three algorithms, its simplicity 
makes it perhaps the best algorithm to use on very dense networks (Exv2). 

Finally, Galil and A. Naamad [12] have recently developed a 
modification to the original Dinic algorithm which results in an algorithm 
with time complexity O({V|-JE|-log7|V|). Their modification of Dinic’s 
algorithm is similar to Galil’s modification of Cherkasky’s algorithm. 
Once again a technique is developed for retaining useful information about 
the structure of the network. 

The remainder of this chapter will be devoted to a closer 


examination of each of the maximum network flow algorithms. 


2.2 Ford and Fulkerson 

Ford and Fulkerson (9) developed the first maximum network flow 
algorithm, known as the labeling algorithm, in 1956. Their algorithm simply 
augments, along some s-t flow augmenting path, the existing tegal flow 
function on a natiisrik. This augmentation is then repeated until there no 
longer exist any s-t flow augmenting paths on the network. In practice the 
zero flow function (i.e. fle)=8 for all edges e) is used as the initial 


existing legal flow function. 
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The labeling algorithm is composed of two basic routines which are 
iterated until a maximum flow function is computed. The first routine 
essentially searches in a systematic way for an s-t flow augmenting path on 
the network. The second routine then augments the existing legal flow 
function along this flow augmenting path. The actual flow augmentation is 
per formed exactly as outlined in the proof of Lemma 1.1. The following 


explanation of the labeling algorithm is taken from (13): 


Step 1. Labeling Process. 

Every vertex is always in one of three states, labeled and scanned, 
labeled and unscanned, or unlabeled. A vertex is labeled and scanned if it 
has a label and we have inspected all vertices adjacent to it. A vertex is 
labeled and unscanned if it has a label but not all vertices adjacent to it 
have been inspected. A vertex is unlabeled if it has no fabel. 

Initially, all vertices are unlabeled. A label for a vertex vj 
always has two parts. The first part is the index of a vertex v;, which 
indicates that we can send flow from vi to Vie and the second part is a 
number which indicates the maximum amount of flow we can send from the 
source to vj without violating the capacity constraints. We first assign 
the label [s*,e(s)=0] to the source, v,. The first label simply says that 
we can send flow from the source to itself; the number » indicates that 
there is no upper: bound on how much can be sent. The source is now labeled 
and unscanned and all other vertices are unlabeled. In general, select a 
ver tex Yj which is labeled and unscanned. Assume vj has a label of the 


form Cit,e(j))] or Ci7,ej)]. For all adjacent vertices vy, Which are 


unlabeled, adjacent to Yj via an edge directed from vj to Vue and for which 
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the edge e=(vi,vy) is useful from vj tow (i.e. flvi vd<clvi syd), assign 
the label [jt,e(k)] to Vv.» Wheres 
e(k) = minle(j), clvi,yd-Flvi vid. 
For ail adjacent vertices YK uhich are unlabeled, adjacent to Vj via an 
edge directed from vy, to Mis and for which the edge e=(vi vj) is useful 
from v5 to vy, (i.e. flvysvj)>8), assign the label! [j~,e(k)] to Ve» where: 
e(k) = minle(j), f(vyyev i). 

The + and the - signs in the labels indicate whether the 
corresponding edges appear as forward or reverse edges in the s-t flow 
augmenting path. Now all the vertices adjacent to Yj have labels; vi is 
considered to be labeled and scanned and may be disregarded during the rest 
of this step. (If one inspects al! the vertices adjacent to Vj and cannot 


J 
scanned vertex.) All the vertices Vv, are now labeled and unscanned. 


label all these vertices, then v. is also considered to be a labeled and 


Continue to assign labels to vertices adjacent to labeled and. 
unscanned vertices until either the sink is labeled or no more labels can 
be assigned and the sink is unlabeled. If the sink cannot be labeled, no 
s-t flow augmenting path exists and, hence, the existing flow function is 
maximum. If the sink is labeled, an s-t flow augmenting path has been 


found and the flow augmentation can be performed using step 2. 


Step 2. Flow Change. 

Assume that the sink is labeled (k*,e(t)]). Let 
flvpavelefly.vylte(t) and turn to vy. If vy, is labeled (jt ,etk)), let 
flvisypdef lv vy d+e(th and turn to vie If vy. is labeled {j°,e(k)], let 


flv, svileflyy.vi)-e(t) and turn to Vie Continue until the source is 
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reached. Erase the labels on al! the vertices and go back to step 1. 


When the labeling algorithm terminates, the set of all labeled 
vertices clearly forms a cut in the network. Further, applying Lemma 1.2 
reveals that the capacity of this cut must be equal to the value of the 
existing legal flow function on the network. Thus by Corollary 1.1, the 
existing legal flow function must be maximum. It now remains to be shown 
that the labeling algorithm will always terminate within a finite number of 
steps. 

Let us first consider the case in which the network is defined by 
an integral valued capacity function. By the same argument as that used in 
the proof of Theorem 1.2, it can easily be seen that the labeling algorithm 
will terminate after at most vif...) iterations, where v(f,,.) is the 
finite value of a maximum flow function on the network. Further, each 
iteration will require at most O(|E]) operations since each edge is 
examined at most twice in the labeling procedure and at most once in the 
augmenting procedure. Thus we have that the labeling algorithm will 
correctly compute, in time OCjEl vif), a maximum flow function on any 
network defined by an integral valued capacity function. It can also be 
shown, however, that there actually exist networks which force the labeling 
algorithm to perform Vit oy! iterations. Consider for example the network 


aT in Figure 2.1.- 
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Network rR 


Figure 2.1 


If the labeling dicoriaki:.Begiantag with the zero flow function on nm. 
augments only along the paths s(s,ala(la,cic(c,flflf,hih(h,t)t and 
s{s,b)b(b,did(d, fi ffc, fice(c,elele,gigl(g,t)t in alternating order, then 2-N 
= VEf os) flow augmentations will be required. Thus the algorithm wil! 
iterate Vito) times. Notice that the inefficiency in this example is 
based on the fact that the labeling algorithm permits the augmentation, at 
each iteration, along any one of several existing s-t flow augmenting 
paths. 

We shal! now consider the case in which the network is not defined 
by an integral valued capacity function. Ford and Fulkerson [3} were able 
to demonstrate the somewhat surprising result that their labeling algorithm 
might fail to terminate if the network was composed of irrational edge 
capacities. This result was based on interpreting the labeling process 
broadly enough to permit the selection of any s-t flow augmenting path at 
each iteration of the computation. Thus the labeling algorithm essentially 
has an unbounded worst case running time. We remark, however, that 
computers only deal with rational numbers and thus in practice we could 


expect the labeling algorithm to halt and yield a correct answer. In fact, 
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despite its weaknesses, the labeling algorithm was successfully used for 


almost 15 years. 


2.3 Edmonds and Karp 

In 1969 Edmonds and Karp [7] showed how the labeling algorithm 
could be modified to obtain a bounded worst case running time. In light of 
our previous remarks, it should not be surprising to learn that their 
modification was essential ly an ordering on the selection of the s-t flow 
augmenting paths. Edmonds and Karp suggested augmenting along the shortest 
s-t flow augmenting path (i.e. an s-t flow augmenting path containing a 
minimum number of edges) at each iteration. This can be easily 
accomplished by modifying the labeling process so that the vertices are 
scanned in the same order in which they receive labels (i.e. by imposing a 
Breadth First Search on the labeling process). The remainder of the 
labeling algorithm is unchanged. The running time bound on Edmonds and 
Karp’s “first labeled, first scanned" modification of the labeling 
algorithm is derived from the following results [15]. 

Consider any network I? upon which there is defined a legal flow 
function f. Let eh) denote the minimum number of edges in an s-u flow 
augmenting path after k augmentations of f. Similarly, let ik) denote the 
minimum number of edges in a u-t flow augmenting path after k augmentations 


of f. 


Lemma 2.1: 


If each flow augmentation of f is made along an s-t augmenting path 
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with’ a minimum number of edges, then: 


a ike (k) 


20) 
and 


Ak) 


(k+1) 
Ty 27) 


for all u,k. 


Proof: (From {[15]) 


(k+1) {k) 


Assume that o <Oy for some u,k. Moreover, let: 
k 
(2.1) gs min fo hth? | g (htt? < a! y, 
Clearly eee 2 1 (only ge ikeh = 8), and there must be some final edge 


(u,v) or (v,u) in a shortest s-u flow augmenting path after the (k+1) 8t 


augmentation of f. Suppose this edge is (v,u), a forward edge, with 


f(v,ul<c(v,u) (the proof is similar for (u,v)). Then @.kth eee) +1 


and by (2.1), 


(k+1) (k) 
(2.2) v, ae saa 1. 


Further, it must have been that f(v,ul=c(v,u) after the th augmentation of 


(k) gg (k) (k+l) 
sz v 


f; otherwise 7, +18), » contrary to the assumption. But if 


f(v,u)=clv,u) after the «th augmentation of f and flv,u)<c(v,u) after the 
(k+1) St augmentation of f, it follows that (v,u) was a reverse edge in the 
(k+1)8t s-t flow augmenting path along which f was augmented. Since that 


path contained a minimum number of edges, 


(k) (k) 
Co aoe Se 1. 


Combining this with (2.2), however, we obtain: 


o'K)42 ¢ g (kel) 


contrary to our assumption. The assumption that gikes} < gin) is 
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therefore false. 


The proof that piste! 2 rik) parallels the above. 


Theorem 2.1: 
If each flow augmentation of f is made along an s-t augmenting path 
With a minimum number of edges, then a maximum flow function is obtained 


after no more than |Vj[-|E[/2 augmentations of f. 


Proof: (From [15]) 

Each time an augmentation of f is made, at least one edge in the s- 
t augmenting path is "critical" in the sense that it limits the amount of 
augmentation. The flow through such an edge (u,v) is either increased to 
capacity or decreased to zero. Suppose (u,v) is a critical edge in the 


(k+1)8* 5-¢ augmenting path. The number of edges in the augmenting path is 


(k) (kK) _ _(k) ({k) 

Co a a SO ea 
The next time edge (u,v) appears in an s-t augmenting path, say the 
(24+1)5t, it will be with the opposite orientation. That is, if it was a 


jst 


forward edge in the (k+1) 8¢, it is a reverse edge in the (f+1 » and vice 


versa. If (u,v) was a forward edge in the (k+1)5* s-t augmenting path (the 


proof is simitar for a reverse edge), then: 


0 fk) = atk) +1 
and 


(2) 
ce ty 1. 


(2) (k) 


(2) ({k) 
By Lemma 2.1, however, ae a “uo = te 


and fT Thus we have that: 
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(2) (2) {k) {k) 
Co ee ea Oe ES 2. 


It follows that each succeeding s-t augmenting path in which (u,v) is a 
critical edge is at least tuo edges longer than the preceding one. 

No flow augmenting path may contain more than |Vj-1 edges. 
Therefore, no edge may be a critical edge more than |V|/2 times. But each 
s-t augmenting path along which f is augmented contains a critical edge. 
Therefore there can be no more than |VJ-|E|/2 successive s-t flow 


augmenting paths and this completes the proof. 


Since Edmonds and Karp’s algorithm differs from the original 
labeling algorithm onty in the order in which the unscanned vertices are 
scanned, it follows that their algorithm will yield a maximum flow function 
if it halts and that their algorithm will require only O(|E]) operations . 
per iteration. By Theorem 2.1, however, Edmonds and Karp’s algorithm is 
guaranteed to halt after at most |V|-|E]/2 iterations. Further, this 
result is independent of the capacity function; holding for trrational 
valued capacity functions as well as integral valued capacity functions. 
Thus Edmonds and Karp’s algorithm will correctly compute, in time 


OUVI-1EI*), a maximum flow function on any network given as input. 


2.4 Dinic 
In 1978, E.A. Dinic [5S] developed a maximum network flow algorithm 
with a bounded OUVI2-JE]) worst case running time. Dinic’s aigorithm, 


like Edmonds and Karp’s algorithm, is based on successive flow 
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augmentations along s-t flow augmenting paths of minimum length. Dinic, 
however, noticed that a single breadth first search of a flow network could 
be used to isolate all minimum length s-t flow augmenting paths. Based on 
this observation, Dinic’s max-flow algorithm is much more efficient than 
the “first labeled, first scanned" algorithm of Edmonds and Karp. In order 
to present Dinic’s algorithm, we must first introduce the notion of a 


layered network. 


Definition 2.1: 

A layered netuwork LN is a network whose vertex set is partitioned 
into disjoint subsets Vg,V,,...,Vg such that Vp = {s} and Vy = {t}. We say 
that V, is the it layer in LN (for @sis&) and that &£ is the length of the 
layered network LN. Each edge e=(u,v) in LN has the property that if ueV; 
then veVi ay (i.e. every sane in LN is directed from one layer to the 


next). 


Example 2.1: (Layered Network LNg) 


Vg={s}. V,={a,b}, Vor{c, d}, V=(t} 


E={(s,a), (s,b), (a,c), (a,d), 
(b,d), (c, t), (d, t)} 


cle)=1 [for alli eck] 


Each iteration of Dinic’s algorithm is called a phase and each 
phase is divided into two procedures. The first procedure generates a 


Jayered netuork from the original input network, in such a way that the 
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layered network isolates all existing minimum length s-t flow augmenting 
paths. The second procedure then uses this layered network to successively 
augment the existing legal flow function along these minimum fength s-t 
flow augmenting paths. We will first describe how Dinic’s algorithm 
performs each of these procedures and then we will state the entire 
algorithm. 

The first task performed during each phase of Dinic’s algorithm is 
the construction of the layered network. The layers composing this network 
are created from a breadth first search of the input network. This breadth 
first search begins at the source and traverses only forward directed edges 
With flow less than capacity or backward directed edges with flow greater 
than zero (i.e. only "useful" edges}. The search terminates when either 
the sink is reached or no nen vertices can be visited and the sink has not 
been reached. When the sink cannot be reached, however, the existing legal 
flow function is a maximum flow function and Dinic’s algorithm halts. 
Notice that performing the search in this way assures that there exists a 
flow augmenting path from the source to each vertex visited. The layers, 
Vie are finally formed by partitioning the vertices visited according to 
the length of their path of discovery from the source. Every vertex vev; 
will then have the property that the shortest s-v flow augmenting path, 
with respect to the existing legal flow function on the input network, is 
of length i. Thereeune; the length of ‘he layered network constructed will 
be equal to the length of the shortest existing s-t flow augmenting path. 


The following procedure formalizes this construction [8]. 
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procedure LN(M, f): 


begin 
Vg: ={s}: 
e i: =@; 


» while Cis=i+l) > 8 do 


‘. begin 
‘ .  Tr=fveVv | vev; for j<i and there exists a 


useful edge from a vertex in V;_} to v}; 
if Te@ then halt (the existing f is maximum) 


else if teT then begin 


‘ Ls =i; 
‘ ‘ ~ — Varmft}s 
. ‘ + ited 
; : end 
. . else Vit=T 
end 
end. 


Once the layers have been generated, the layered network’s edge set 
and capacity function are constructed. The edge set is simply constructed 
from all edges in the original input network which are useful from a vertex 
in V; to a verhes in Vii, (for all @<is£). Every edge in the layered 
network, however, is directed from the jth layer to the (141) 8¢ layer, 


regardless of its orientation in the input network. Any edge in the 


layered network whose orientation is different in the input network is said 
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to be a reverse edge in the layered network. Similarly, any edge in the 
layered network whose orientation is the same in the input network is said 
to be a forward edge in the layered network. The capacity function C is 
then created from the excess capacity along each edge in the layered 
network as follows: 
Clelecle)-fle), for all forward edges in the layered network 
and 
cledefle), for all reverse edges in the layered network, 

where f is the existing legal flow function and c is the capacity function 
on the input network. [It should now be clear that every s-t chain in the 
fayered network corresponds to an existing minimum length s-t flow 
augmenting path in the input network. Further, every minimum length s-t 
flow augmenting path in the input network is represented by an s-t chain in 
the layered network. Thus the layered network constructed during each 
phase of Dinic’s algorithm essentially isolates all existing minimum length 
s-t flow augmenting paths. 

The second task performed during each phase of Dinic’s algorithm is 
the flow augmentation. Starting at the sink in the newly constructed 
layered network, the algorithm follows edges backward to the source to find 
an s-t flow augmenting path. (Recall! that every s-t chain in the layered 
network corresponds to an existing s-t flow augmenting path in the input 
network.) The existing legal flow function is then augmented along this 
path as outlined in the proof of Lemma 1.1. Notice that the excess capacity 
along each edge in the flow augmenting path is given by the layered network 


capacity function, C. After augmenting along this path, the algorithm 
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adjusts the capacity function c to reflect the new excess capacity along 
each edge in the path. It also deletes from the layered network all edges 
in the path whose excess capacity drops to zero. Finally, it deletes from 
the layered network all vertices and their incident edges which are no 
longer reachable from the source or the sink. This is accomplished by 
deleting all vertices which either have no incoming edges or have no 
outgoing edges, and continuing to delete such vertices until every vertex 
left in the layered network has at least one incoming and one outgoing 
edge. After all necessary deletions have been performed, the algorithm 
searches for another flow augmenting path and the augmentation process is 
repeated. This continues until there no longer exist any s-t chains in the 
layered netuork. 


Dinic’s complete max-flow algorithm is: (From (28)) 


procedure DINIC(/Q): 
begin 
. initialize existing legal flow function f, on input network 7M, to @; 


» while “true” do 


. » construct layered netuork, LN; 
. - for each vertex v in LN do 

- : - calculate indegree (v); 
‘ ‘ - calcutate outdegree (v); 


P , - if Cindegree (v)=8) or foutdegree (v)=8) then 


, . ; add v to nullist 
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. . end; 
. - while t is a vertex in LN do 
‘é - begin 


: ‘ . trace back from t to s to find an augmenting path; 


. : + augment f along this path; 

. . - update c with new excess capacity along each edge in path; 
. - - delete from LN all edges along path which now have zero 
: é F excess capacity (i.e. cle) =8), updating indegrees, 

. . ; outdegrees, and nullist; 

: . » while some vertex v is on nullist do 

‘ : . delete v and incident edges from LN and from nullist, 
- . . updating indegrees, outdegrees, and nullJist 

. : end 

/ end 

end. 


Recall that Dinic’s algorithm terminates when the breadth first 
search performed in constructing each layered network fails to reach the 
sink. When this occurs, however, the set of all visited vertices clearly 
forms a cut in the input network. Further, applying Lemma 1.2 reveals that 
the capacity of this cut must be equal to the value of the existing !egal 
flow function. Thus by Corollary 1.1, the existing legal flow function 
must be maximum. It now remains to be shown that Dinic’s algorithm will 
always terminate within time OUVI2-JE}). This is proven as a consequence 


of the following lemma [8] which shons that the number of phases is bounded 
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by |VI. 
Let a denote the length of the |ayered network constructed during 


;th 


the phase of Dinic’s algorithm. 


Lemma 2.2: 


If the (k+1)%* phase of Dinic’s algorithm is not the last, then 


esl ig o,. 


Proof: {From [8)) 
Consider any s-t chain in the layered network constructed during 


the (k+1)8t phase of Dinic’s algorithm: 


a heat 
——— 
> ss a as PE t. 


First, let us assume that all the vertices in this chain appear in 
the kth layered network. Let vi be the jth layer of the wth layered 
network. We claim that if VaeVp then a2b. This is proven by induction on 
a. For a=8, (vg=s) the claim is obviously true. Now assume v,,,€V_.- If 


csb+l the inductive step is trivial. If, however, c>b+l then the edge C541 


«th 


was not used in the «th phase since it was not even in the |ayered 


network, in which only edges between adjacent layers appear. But if C54) 


was not used in- the kth phase and is useful from v, to in the 


a Vatl 


beginning of the (k+1) 8t phase, then it was useful from v. to v in the 


a atl 


beginning of the Kth phase. Thus, Vag] cannot belong to Ve (by procedure 
LN). Now, in particular, sa TA 1 and ae Therefore, Ly 4122, 
+ 


Further, equatity cannot hold because then the entire s-t chain would have 
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been’ in the th layered network, and if all its edges are still useful at 
the beginning of the (k+1)8¢ phase then we have a contradiction to the 
th 


termination of the phase. 


If not all the vertices in the s-t chain appear in the th layered 
fa+l 


network then let Va Yagy bE the first edge such that for some b, 


Val, but v is not in the kth jayered network. Thus, G54) as not used 


atl 
in the «th phase. Since it is useful in the beginning of the (k+1) 8¢ 


th phase. Thus 


phase, however, it was also useful in the beginning of the 
the only possible reason for v,,,; not to belong to V,,, is that b+i=2,. 


Further, by the argument of the previous paragraph azb. Therefore, at+l2£, 


and so Dy 41> 4 


Corollary 2.1: 


Given any network f as input, the number of phases performed by 


Dinic’s algorithm must be less than or equal to |V|. 


Proof: 
Any layered network constructed by Dinic’s algorithm must contain 
no more than |V| layers. Thus by Lemma 2.2 there can be at most |V] 


phases. 


0 


We now notice that the time required during each phase of Dinic’s 
algorithm, to construct the layered network and initialize the indegrees, 
outdegrees, and nullist is bounded by O(jE|). Further, each flow 


augmentation requires time O({V|) and there can be at most O(|Ej) such 
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augmentations since each augmentation causes the deletion of at least one 
edge from the layered network. Finally, the total time required during 
each phase of Dinic’s algorithm to davete edges, delete vertices, and 
update indegrees, outdegrees, and nullist is bounded by O(jEJ). This 
results from the fact that each edge and each vertex can be deleted from 
the layered network at most once. It should now be clear that each phase 
(iteration) of Dinic’s algorithm has time complexity O(|V]-|E]). Thus by 
Corollary 2.1, we have that Dinic’s algorithm will always terminate within 
time O(V|2-|E)). Therefore, Dinic’s algorithm will correctly compute, in 


time OUVI2-JE]), a maximum flow function on any network given as input. 


2.5 Karzanov 

A.V. Karzanov [8,14] modified Dinic’s algorithm in 1973 to obtain 
an ojvys) maximum network flow algorithm. Karzanov noticed that the 
layered network constructed during each phase of Dinic’s algorithm could be 
used to simultaneously augment along all existing minimum length s-t flow 
augmenting paths. He then showed how this simultaneous augmentation could 
be performed in time OcVv|2). Karzanov’s algorithm is the result of 
replacing Dinic’s OC{V|-JE]) successive flow augmentation procedure with 
this new Ov?) simultaneous flow augmentation procedure. 

Karzanov’s results are based on the notion of a maximal flow 
function. A santa flow function t, on a network It, is defined to be any 
lega! flow function on Mf which has the property that every s-t chain in It 
contains at least one saturated edge (i.e. at least one edge e such that 


f(el=c(e)). From the fallowing example it can be seen that a maximal flow 
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function on a network need not be a maximum flow function on that network. 


Example 2.2: 


The following maximal flow function tg is defined on the layered 
network LNg from Example 2.1: 


fg(s,a)ql Fg(s,b)=8 


fglacl-@  fgla,dd=l 


fglb,d)=8 fg(c, t)=8 


fg(d, t) =] 


v (Fg) el is not maximum. 


Karzanov noticed that the flow augmentation required during each 
phase of Dinic’s algorithm could be achieved by simply augmenting the 
existing !egal flow function f with any maximal flow function # on the 
current layered network. Once such a maximal flow function had been 
| computed, the flow modification could be performed as fol lows: 

t’(elefle)+f le), for all forward edges in the layered network 

and 

f’(elefle)-fle), for all reverse edges in the layered network. 
It can easily be seen that performing the augmentation in this manner is 
essentially the same as simultaneously augmenting along all existing 
minimum length s-t flow augmenting paths, where each individual 
augmentation is performed exactly as outlined in the proof of Lemma 1.1. 
Karzanov used this clever augmentation technique to develop the following 


modification of Dinic’s algorithm. 
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procedure KARZANOV (72): 
begin 
- initialize existing legal flow function f, on input network I, to 8; 


-» while "true" do 


. begin 

. - construct layered network, LN; 

. » compute maximal flow function # on LN; 

. - for each edge e in Ido 

. . if e is a forward edge in LN then fle)r=fled+fle) 

. . else if e is a reverse edge in LN then fle) :=f le) -# le) 
- end 

end. 

Notice that Karzanov’s algorithm, like Dinic’s algorithm, 


terminates when the breadth first search performed in constructing each 
layered network fails to reach the sink. By the same argument as that used 
for Dinic’s algorithm, however, the existing legal flow function must be 
maximum when this occurs. Therefore if Karzanov’s algorithm halts, then 
the existing legal flow function must be a maximum flow function on the 
input network. Next we notice that the proof of Lemma 2.2 is valid for 
Karzanov’s algorithm as well as Dinic’s algorithm... Thus given any network 
Mas input, the number of phases performed by Karzanov’s algorithm must be 
less than or equa! to |V|. It should also be clear, however, that the time 
required to construct each layered network and perform each flow 
modification is bounded by O(JE|). Therefore if we let ¢ denote the time 


required to compute each maximal flow function, then Karzanov’s algorithm 
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is guaranteed to halt within time O(|V{-{E] + [V]-2). We will now show hon 
Karzanov’s algorithm computes each maximal flow function in time o1v}?) to 
yield an O({V|3) max-flow algorithm. 

Karzanov’s algorithm computes each maximal flow function by 
successively improving an existing illegal flow function, called a preflow 
function, on the current layered network LN [8,18,14). For each vertex v 
in LN, let In’(v) denote the set of all edges incoming to v in LN and let 
Out’(v) denote the set of all edges outgoing from v in LN. A preflow 
_ function cs on the layered network LN, sence Stes with each edge e in LN a 


real number fle) satisfying the conditions: 


C3) Q<fle)scle), for each edge e in LN 
C4) > fle) 2 >> fle), for each vertex vz#s,t in LN. 
eeln’ (v) eeOut’ {v) 
Thus a preflow function is simply a flow function which satisfies the 
capacity constraint but not necessarily the conservation constraint of a 
legal flow function. For each vertex v#s,t in LN, we define excess(v) to 
be the excess flow entering v: 


excess{v) = Dd fle) - D fled. 
eeln’ (v) eeOut’ (v) 


If excess(v)>8 then v is said to be unbalanced; otherwise v is said to be 

balanced. The source and the sink are always considered balanced. 
Throughout the execution of Karzanov’s maximal flow procedure, 

every edge in LN is declared either open or closed. Initially all edges 


are declared open. As the algorithm proceeds, however, some of the edges 
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Will’ be declared closed. Once an edge is declared closed, the flow through 
it Hill onctts unchanged to the end of .the procedure. 

Karzanov’s maximal flow procedure alternates between pushing 
additional flow from unbalanced vertices and balancing the unbalanced 
vertices that are generated during these pushes. The pushing of flow is 
achieved through repeated calls to a procedure PUSH(i), with increasing i 
(18]. The procedure PUSH(i) considers in turn each unbalanced vertex in 
layer V;, attempting to push flow from it to vertices in layer Viate For 
each unbalanced vertex veY;, the procedure considers in turn each open 
edge in Out’(v) and sends through it the maximum possible amount of flow. 
(The two constraints that exist are the current excess of v and the amount 
of flow needed to saturate the edge.) The push from v ends when either v 
becomes balanced or every edge in Out’({v) becomes either saturated (i.e. 
#(electel) or closed. For each vertex u in LN there is a stack (push-down 
stored on which the history of additions of incoming flow into u is 
recorded. When the flow in an edge e=(v,u) is incremented by an amount 6, 
the pair (v,6) is added to the top of the stack for vertex u. The procedure 
PUSH(i) is said to be successful if flow is pushed to layer Vigte lf 
PUSH(i) is successful then PUSH(i+1) is called, and so on. 

A procedure BALANCE(i) is the tool through which unbalanced 
vertices become balanced. The procedure BALANCE(i) uses the stacks to 
shift back flow fen vertices in layer V; to vertices in layer Vi_qe It 
balances im turn each of the unbalanced vertices veV; by canceling the 
most recent additions of flow into v. Clearly the last canceled addition of 


flow into v may only be partial. After each unbalanced vertex v is 
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balanced, all the edges in In’(v) are declared closed. The procedure 


BALANCE (i) is always followed by a call to PUSH(i-1). 


Karzanov’s complete maximal flow procedure is: (From (18]) 


procedure MAXIMAL (LN): 


PLOOP: 


initialize existing preflow function t on layered network LN, to 8; 
empty the stacks of all vertices in LN; 

i:=@; 

PUSH (i); 


. while the previous push was successful and i+l<2 do 


begin 
: ‘ is=i4l; 
: .  PUSH(i) 
: end; 


if there exist unbalanced vertices in LN then 
begin 
is=number of highest layer V ; (8<j<2) containing 


unbalanced vertices; 


BALANCE (i); 
. is=i-1; 
: - -goto PLOOP 
end; 


for each edge e in LN do fle): =f (e) 
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It should now be clear that every vertex in the layered network LN 


Will be balanced when Karzanov’s maximal flow procedure halts. Thus the 
final existing preflow function f will in fact be a legal flow function on 
LN. In order to show that it will also be a maximal flow function on LN, 


we must introduce the notion of a blocked vertex. A vertex v in LN is said 
to be blocked with respect to the existing preflow function if every v-t 
chain in LN contains at least one saturated edge. Notice that the source s 
becomes blocked after the first execution of PUSH(®8), since every edge in 


Out’ (s) becomes saturated. 


Lemma 2.3: 
If a vertex in LN becomes blocked at some point in the execution of 
Karzanov’s maximal flow procedure, then it remains blocked to the end of 


the procedure. (A proof of this lemma appears in [8].) 


Lemma 2.4; 
Every vertex in LN is balanced at most once throughout the 


execution of Karzanov’s procedure. (A proof of this lemma appears in [8].) 


We can now see that the final existing preflow function on LN will 
be a legal flow function which has the property that every s-t chain in LN 
contains at idaat one saturated edge. Therefore when Karzanov’s maximal 
flow procedure halts, the existing preflon function ¢ will in fact be a 
maximal flow function on the layered network LN. It now remains to be 


shown that the procedure will always halt within oCvy2) steps. 
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The total number of steps performed by Karzanov’s maximal flow 
procedure is clearly bounded by the total number of flow additions and flow 
reductions performed. The number of flow reductions performed, however, is 
bounded by the number of flown additions performed since each vertex is 
balanced at most once and the history of flow additions in the stacks is 
used to perform the flow reductions. Thus it suffices to show that the 
number of flow additions performed by Karzanov’s procedure is bounded by 
O(yVy2). We first notice that there can be at most one saturating flow 
addition per edge in the layered network. Since the number of edges in the 
layered network is bounded by the number of edges in the original input 
network, however, there can be at most O({E[) saturating flow additions. 
Next we notice that there can be at most one non-saturating flow addition 
per vertex in the layered network, between any two successive calls to 
BALANCE (i). (When flow is pushed from a vertex v in layer Vie only the 
last edge considered in Out’{v) does not necessarily become saturated.) 
Since the number of vertices in the layered network is bounded by the 
number of vertices in the original input network, however, there can be at 
most O(|V]) non-saturating flow additions between any two successive calls 
to BALANCE(i). From Lemma 2.4, however, there can be at most O(|V|) calls 
to BALANCE (i). Thus the total number of non-saturating flow additions is 
bounded by o(yv}2) and hence the total number of flow additions is bounded 
by O(|V{2). | 

It should now be clear that Karzanov’s maximum network flow 
algorithm will correctly compute, in time O({V|-JE] + IV}-1V 12) = o1v3), a 


maximum flow function on any network given as input. 


49 


2.6 Further Improvements 

We saw in the last section Hakan O(t) algorithm for computing a 
maximal flow function on a layered network could be used to develop an 
OCUVI-TE] + |V]-t)) maximum network flow algorithm. (Simply modi fy 
Karzanov’s max-flow algorithm by replacing his "push and balance" procedure ~ 
with the new maximal flow procedure.) Notice, havevers that any max-flow 
algorithm developed through this technique can be no faster than 
OC(|VI-[E{). Despite this fact, each of the four most recent maximum 
network flow algorithms have been based on developing new maximal flow 
procedures. We shall now briefly describe each of these new maximal flow 
procedures. 

In 1976 B.V. Cherkasky [18] showed how Karzanov’s oy?) push and 
balance routine could be modified to run in time ocpy| epi’). 
Cherkasky’s procedure partitions the layered network into blocks of 
consecutive layers called superlayers. It then applies Karzanov’s push and 
balance techniques to these supertayers. Within the supertayers, however, 
Dinic’s flow augmentation techniques are used. The result is an 
asymptotically faster but very complex maximal flow procedure. 

A tittle over a year later, Z. Galil (18) improved Cherkasky’s 
routine to obtain a maximal flow procedure with time complexity 

o(|v{2/3.4E 52/3) Galil’s procedure differs from Cherkasky's procedure in 
the techniques used within the superlayers. Galil’s routine maintains a 
special data structure containing information about the current 


"usefulness" of chains within the layered network. This data structure is 
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used to expedite the push of flow through edges within the super layers. 
Like Cherkasky’s routine, Galil’s procedure is very complex and requires a 
great deal of overhead. 

In 1978 a very simple O{v|) maximal flow procedure was developed 
by three Indians named V.M. Malhotra, M. Pramodh Kumar, and S.N. 
Maheshwari [16]. Their procedure is based on successively augmenting an 
existing legal flow function on the layered network. The procedure begins 
by determining the maximum amount of flow that can be pushed through each 
vertex v in the layered network. This value is called the flow potential 
pb,iv) of the vertex v. Each flow augmentation is then performed in three 
steps. First a vertex with minimum non-zero flow potential over all 
vertices in the layered network is selected as the reference vertex, r. 
Next, p,(r) units of flow are pushed from r to t and from s to r. The 
pushing is performed essentially as outlined in Karzanov’s algorithm. 
Finally, the procedure updates the flow potential of each vertex through 
which flow has been pushed, closing all edges which become either saturated 
or unreachable from s or t. Although this procedure is not an asymptotic 
improvement over Karzanov’s push and balance routine, it is extremely 
simple and results in perhaps the best max-flow algorithm for use on dense © 
netunorks (JE, =1V[2). 

Finally, Galil and A. Naamad [12] have recently developed an 
OtiELAcgeWilvdaxiaal Flea lproceduces - This: proceduke 1s. eimiiak te 
Galil’s o(jvye/3.16 12/3) maximal flow procedure in that a data structure is 
maintained for processing chains within the layered network. The new 


algorithm, however, does not partition the layered network into 


51 


super layers. 

This concludes our discussion of algorithms for computing a maximum 
flow function on a network. We will now turn our attention to the question 
of lower bounds on the soncutat ional complexity of the maximum network flow 


problem. 
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CHAPTER 3 - LOWER BOUNDS ON MAX-FLOW 


3.1 Introduction 

Chapter 2 dealt with the establishment of upper bounds on the 
computational complexity of the maximum network flow problem. We traced 
the development of max-flow algorithms from the original labeling algorithm 
of Ford and Fulkerson [9], through the recent OV] JE} -log@{V}) algorithm 
of Galil and Naamad [12]. The very fact that the search for new max-flow 
algorithms has been so fruitful leads us to now ask the question, “Can we 
do better?". Can we develop a max-flow algorithm which is asymptotically 
faster than OUIV|-[E|-log@}V¥])? Galil has shown [11,12] that any algorithm 
which uses Dinic’s technique of dividing the problem into phases (as do all 
the known algorithms developed since Dinic’s algorithm) must have time 
complexity at least O({V|-jE|). Further, he has conjectured an N(|V|-JE]) 
lower bound on the computational complexity of the maximum network flow 
problem. At this time, however, there is no known non-linear lower bound 
on max-flow. 

The determination of lower bounds on the computational complexity 
of a problem is generally much more difficult than the establishment of 
upper bounds on the problem. In the fatter case we can simply demonstrate 
an algorithm for solving the problem within the specified running time. In 
the former case, however, we must prove that any algorithm for solving the 
problem must require at least the specified running time, regardless of how 
clever the algorithm. In order to somewhat simplify the lower bound 


problem, many authors have chosen to work with restricted models of 
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computation. One such mode! which has received considerable attention in 
the recent literature is the linear decision tree model (6,18,21,22]. This 
mode! tends to underestimate total time complexity but nevertheless enables 
us to study non-trivial lower bounds. In this chapter we shall investigate 
one particular approach to establishing non-linear lower bounds on the 
computational complexity of the maximum network flow problem relative to 
the linear decision tree mode! of computation. The technique we shall deal 
With is the polyhedral technique developed by A.C. Yao, O.M. Avis and R.L. 


Rivest [21,22]. 


3.2 The Model of Computation 

The linear decision tree model of computation is based on the 
notion of a linear decision tree algorithm. A linear decision tree 
algorithm, operating on input (xp eee KA) is simply a finite ternary tree 
with each internal node representing a test of the form "Z @,-x,2z" and 
each leaf containing a possible output. Given any input, the algorithm 
begins at the root and proceeds by moving down the tree until a leaf is 
reached. At each internal node the algorithm performs the specified test 
and then branches according to the result of this test (<,=*,or >). Once a 
leaf is reached, the information contained in that leaf its output as the 
result of the computation and the algorithm halts. The time complexity of 
any such aigor ttbn is simply defined to be the height of the corresponding 
tree. 

The computational complexity of any problem relative to the linear 


decision tree mode! of computation can now be defined as the minimum height 
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over all decision trees which solve the problem. It should be clear that 
the linear decision tree model measures time complexity solely in terms of 
the number of comparisons and branchings required. Thus the model tends to 
underestimate the total time complexity of a problem. Despite this fact, 
the linear decision tree model has proven useful in establishing several 


non-trivial lower bounds. 


3.3 Polyhedral Decision Problems 
Let P = fer” | 2,608 for each i=1,2,...,m} be a set of points in 


R", where Relxpreeer Xp), m is an integer and 


n 
j= 


for real numbers hij: The set P is said to be a polyhedron in mR", If d 
is the dimension of the smallest subspace of R™ containing P, then P is 
also said to be a polyhedron of dimension d. (Notice that we are 
restricting our attention to homogeneous polyhedra, i.e. cones.) On each 
subset H of the set {1,2,...,m}, we define the set of points Fy(P)cP as 
fol lous: 
Fy(P) = {xeR” | 2: (2)<8 for each ieH, 
2, (28 for each i ¢H}. 

The set Fi,(P) is called a face of the polyhedron P. If s is the dimension 
of the smallest subspace of R" containing FitP), then Fy(P}) is said to be 
a face of dimension s. (The empty face has dimension -1 by convention.) We 
shall tet F.(P) denote the set of all faces of dimension s of P. Notice 


that every point in P lies on some face of P and that the intersection of 
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any two faces of P is empty. 

The polyhedral decision problem B(P) can now be defined as the 
problem of determining whether an input point ReR” fies in the polyhedron 
P (i.e. Given any input %eR™, is XeP?). A linear decision tree algorithm 
for solving this problem will be a decision tree which contains a "yes" or 


no" decision at every leaf. The computational complexity of B(P) relative 
to the linear decision tree model of computation will be denoted by L(P). 
In 1977 A.C. Yao, O.M. Avis and A.L. Rivest [21] proved the 


following fundamental theorem relating the complexity of the polyhedral 


decision problem B(P) to the facial character of the polyhedron P. 


Theorem 3.1: 
Let P = {% | 2,(%)<@ for each i=1,2,...,m} be a polyhedron in R". 
Then for each s, 


L(P) > 1/2 log IF, (P|. 


Proof: 


The proof of Theorem 3.1 can be found in [21]. 


This theorem states that Q(log Fo) linear comparisons are necessary to 
determine if a point lies in a polyhedron composed of Fs s-dimensional 
faces. As a result, we can determine lower bounds on the computational 
complexity of any polyhedral decision problem relative to the linear 
decision tree model of computation by simply examining the facial structure 


of the polyhedron. 
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3.4 Applications to Max-Flow 

In this section we shal! consider a straightforward application of 
the concepts presented in the previous section to the problem of 
establishing a non-trivial lower bound on max-flow. We shall first 
introduce the class of polyhedral decision problems {B(P,) | n22}, which is . 
very closely related to the maximum network flow problem. We shall then 
formalize this relationship by showing that L(P_)-n2-2 is in fact a lower 
bound on L,, where L, is the linear decision tree complexity of the maximum 
network flow problem for a complete network on n vertices. Thus any lower 
bound on L{P.) will also yield a lower bound on L,. Finally we shall prove 
the following three results concerning the facial structure of the 
polyhedra P, (n22): 


1) There exists a positive constant c such that 
IF, (PDI 2 c(n-2)!, for all n22. 


2) There exists a positive constant c’ such that 
IF, (P)] < c%+(n-2)!, for all n22. 


3) There exists a positive constant c’” such that 
IF PL s 2°77), tor alls and for all n22. 
Based on these results, we can then conclude that Theorem 3.1 cannot be 
directly applied to the class of problems {B(P,) | n22} to obtain a non- 
trivial tower bound on max-flow. 

In order to formally define the class of polyhedral decision 
problems {B(P.) | 22}, we must introduce some new notation. Let G, denote 
the complete directed graph on n22 vertices in which the vertex v, is 
specified as s and the vertex v_ is specified as t. We can represent the 


n 


sets of vertices and edges of G., as follows: 


V = {vy .Vo.Vgee00s¥p} 
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and 
E = fe; ; | 1si, jsn}. 


Notice that |V]=n, [E|=n@ and we are defining Gj to be the edge in CG. 


directed from vertex v; to ver tex vj It should be clear that any capacity 


function defined on the set E will give rise to a complete flow network on 
n vertices. Further, notice that the definition of a cut in a network can 
be directly applied to the graph G,. We now let {Xp Xoe0+-Xoma} represent 


the set of al! cuts in Gh 


2 
For our purposes, a vector Jer" +1 Hill be represented as 


+ ‘ n?4+1 
Ge (Up pee ees U gre ee Unpee s+ *Unn nay)» The polyhedron P, in R can now 


be defined as follows: 
a 2 ne: -2 
P= {GeR" +1 | yj j28 for all 1si,jsn, yn24128, 4.(9)28 for all 1sk<2"%}, 


uhere 


4 = (Z ty; ; | 2) 63%) -Uneat 
The polyhedral decision problem B(P,) is to determine whether an input 


2 
point der" +1 belongs to the polyhedron Pye If we think of the set 


{Urq++++*Unnd a8 defining a capacity function on the graph G, (i.e. 


2 
cle; j)=u;j). then it should be clear that a point Ger” +1 will belong to 


the polyhedron Pa if and only if the following two relations hold: 


1) y; 28, for all isi, jsn (i.e. cle; ey; j defines 
a legal capacity function on G)). 


2) Bsy,2,1<C(X%,;,), where X.., is any minimum 
capacity cut on the network defined by G,, 


and cle; j)=u;j- 


By Theorem 1.3, however, we have that CX, ip! is equal to the value of any 


maximum flow function f,,. on the network defined by G,, and cle; =u; 5. 


x 
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Therefore, it should also be clear that a point jer” + will belong to the 
polyhedron P. if and only if cle; qu; ; defimes a legal capacity function 
on G, and there exists some legal flow function f on the network defined by 
G, and cle; j)=u;j- such that v(fl=u,2,) (i.e. OSu,2,,SVlfma,d)- The 
following lemma relates the linear decision tree complexity of BI(P,) to the 


linear decision tree complexity of the max-flow problem on a complete 


network of n vertices. 


Lemma 3.1: 


L, 2 L(P,)-né-2. 


Before presenting the proof of Lemma 3.1, we should consider the 
structure of any linear decision tree algorithm which computes a maximum 
flow function on a complete network on n vertices. Such an algorithm will 
be a ternary tree operating on input (Up pee Un)» Where the input defines 
a capacity function cle; uj; on the graph G.. Each leaf in the tree Hill 


contain a set of n@ 2 


linear functions (gi; | 1si,jsn} defined on the n 
input variables. For any input, the algorithm will begin at the root and 
proceed by moving down the tree until a leaf is reached. Once a leaf is 


reached, a maximum flow function fmax on the network defined by the graph 


x 
G, and the input (yjyreeesUgn). Will be given by fmax (ej j) = 
Di jY1pee e+ Un) We should note that on every network there will exist a 


maximum flow function that can be completely defined by a set of linear 


combinations of the edge capacities on the network. 
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Proof of Lemma 3.1: 

Let T be any optimal linear decision tree algorithm for computing a 
maximum flow function on a complete network on n vertices. Clearly the 
height of T must be L,. Further, we can obtain a linear decision tree rs 
for the problem B(P,.) by modifying T as follows. Place the root of T below 
@ new sequence of né distinct tests of the form "Is yj j28?" such that the 
root of T is reached if and only if all of these new tests produce a "yes" 
answer. Then replace each leaf in T with a new test of the form "Is 
Up24129?" followed by a new test of the form "Is V (Ff nay) “Yn241 29?" (Yno4) 
being a new input). The new tree T” should be constructed in such a way 
that if any of the newly added tests produce a "no" answer, then a leaf 
containing a "no" decision is reached. Otherwise a leaf containing a "yes" 
decision is reached. Since the value VA Fiax? is simply a sum of the ij 
available at each leaf in T, it should be clear that T” is in fact a linear 
decision tree algorithm for solving the problem B(P,). The height of T’, 


242 and thus we have the relation: 


however, is Lytn 
LIP.) < Lyen?42 


=> L, 2 LIP.)-n@-2. 


The problem of establishing a non-trivial lower bound on the linear 
decision tree ecupientty of the maximum network flow problem has now been 
reduced to the problem of establishing a non-trivial lower bound on L(P,), 
for each n22. By Theorem 3.1, however, we can establish lower bounds on 


L(P.), for each n22, by simply examining the facial structure of the 
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polyhedra Pas The remainder of this section will be devoted to proving 
several lemmas concerning the number of faces composing each of these 
polyhedra. We will essentially show that each polyhedron is composed of 
relatively few faces and thus Theorem 3.1 can be of no use in establishing 


a non-linear lower bound (if one exists) on any L(P,). 


Lemma 3.2: 


There exists a positive constant c such that JF, (PA) I 2 c(n-2)!, 


for all n22. 


Consider any n22 and let p” be an s-t chain in the graph G.. Now 
consider the point Bert! with the properties: 
1) uj j= if ej eR’ for all lsi,jsn 
2) yj j78 if e; #P for all 1si,jsn 
3) Un24i = 
Clearly the point §” belongs to the polyhedron re and so there must exist 
some face of P,, containing GU’ (since every point in a polyhedron lies on 
some face of that polyhedron). Let Fig (P) denote the face of P, which 
27 


contains the point UY’. Therefore, Fy: (P,) is the set of al! points U’’eP 


which satisfy: 


(3.1) yj j=8 iff yj jn for all isi, jsn 
(3.2) u75 4128 
(3.3) O98 iff 0-8, for all 1<k<2-2, 


From (3.1) and (3.2) it should be clear that every point U*“eFy- (P,) will 


be non-zero in exactiy the same co-ordinates. Further, applying this 
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observation to (3.3) and the fact that p’ is an s-t chain in G,, we must 
have that Uj j"Yn24l for each non-zero Uij (1<i,jsn). Thus it must be the 
case that every point in the set Fa (P,) is simply a scalar multiple of the 
point 9”. It should now be clear that the smallest subspace of Reel 
containing the set Fg (P,) has dimension 1. Thus we have that 
Fa-(P)eF, (PL). Finally, notice that each distinct s-t chain in G, will 
give rise to a distinct 1-dimensional face of P.. Therefore, if Z denotes 
the number of distinct s-t chains in the graph G, then we must have: 

(3.4) IFy (P| 2 2. 


Since Gh is a complete graph on n vertices, however, the value Z can be 


expressed by the following formula: 


n-l 
Z= D> (n-2)1/(i-1)! (8! = 1) 
i=l 
n-2 
= (n-2)!-)) 1/1! 
i=6 
> (n-2)! 


Combining this expression with (3.4) we obtain: 


IF, (PV) 2 (n-2)! 


Lemma 3.3: 


There exists a positive constant c” such that JF, (P,)] < c*(n-2)!, 


for all n22. 


Proof: 


Consider any n22 and let Fi (P.) denote the set of all 1-dimensional 
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faces of the polyhedron Pe which arise from distinct s-t chains in the 
graph G., as outlined in the proof of Lemma 3.2. Now let Q denote the set 
of all points U’eP,, which satisfy any one of the following three 
conditions: 

1) 9-0 

2) Un241 78) Y{p_>z8 and Y” has exactly one non-zero co-ordinate 

3) YU’ belongs to a face in the set Fr (P.). 
First notice that each point in Q which satisfies condition (2) belongs to 


2 


one of exactly n°-1 distinct (trivial) 1l-dimensional faces of P.. Further, 


notice that none of these n@ 


-1 faces belongs to the set Fy(P.). Thus each 
point in Q, except the point 8, belongs to one of exactly |F{(P) |4n2-1 
distinct l-dimensional faces of P+ Next observe that every point in the 
polyhedron P, can be expressed as a convex combination of points in Q 
(consider expressing a network as a sum of distinct s-t chains and isolated 
edges). Therefore, every point on a l-dimensional face of P, can be 
expressed as a convex combination of points in Q. Thus we must have that 
the set Q contains at least one point belonging to each 1l-dimensional face 
of P,, and hence: 
IF, (PL)] s IF((P | + n% - 2 


n-2 
= ((n-2)!-D isi!) + n*% - 1 


i= 
< 3-(n-2)! + n% - 1 


< c%+(n-2)! (for some positive constant 
c’, independent of n) 


| 
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Lemma 3.4: 


There exists a positive constant c’” such that |F,(P,)] s 


se, 2 
2'er-n for all s and for all n22. 


Proof: 
2 
Consider any n22 and let each vector eRe” be represented as 
; 2 
= (Wy peers Mage Ayp seer yA ay) A polyhedron es in R°™ can now be defined 
relative to the graph G. as follows: 
2 
P= feRe” | HW; ;28 and Hi; 28 and Wy 2H for all 1si, jsn, 
fe (B28 and Be (A) 8 for all 2sqsn-1, 
B28 for all isks2"-4}, 
where 
47> _— — 
ao = D fa; | ej jeln(v)} - D (8; | e; jedut(v,)} 
and . 
Oe = DY fay, Ley eek} - DL fj; | ej jelntt)}. 
If we think of the set {pyre eee Hand as defining a capacity function on the 
graph CG, (i.e. cle; =n; 5) and the set {yp poeee Aaah as defining a floun 
function on the network defined by G, and cle; dan; | ij 
2 
then it should be clear that a point eRe” will belong to the polyhedron 


(i.e. fle; =u ), 


ps if and only if the following two conditions hold: 


1) cle; ,)=w,, defines a legal capacity function 
“on the graph G,, (i.e. cle)28 for each eeE). 


2) fle; )J=n;, defines a legal flow function on 


the network defined by G, and cle; i) =; |. 


Notice that the constraints Oy” (i) 28 are all redundant since any legal flow 


function f on a network Tt must always satisfy the relation v(f)<C(X,; 7), 
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where X..,. is any minimum capacity cut on Mt. It should now be clear that 


2 
we can determine whether a point eRe” belongs to the polyhedron P by 


242n-4 constraints 


simply testing to see if i satisfies each of the first 3n 
defining Pr. Thus there exists a straightforward linear decision tree 
algorithm of height 3n%42n-4 for solving the problem B(P?) and hence we 
have the relation: 
LIP“) < 3n%42n-4 < Gn?. 

Further, combining this relation with Theorem 3.1 we obtain the result: 

4n@ > LIP’) 2 1/2 log IF,(P4)|, for all s 
(3.5) = 28n" , IF (P|, for all s. 
If we now let F(P?) = U F(Pi) denote the set of all faces of the 


polyhedron P/, then by (3.5) we have: 


2 
FP] < > 289 


-1s8<2n¢ 


= |FIPZ)| < (2n242)-28"" 
(3.6) = |FIP)| < 2fe°n") 
(for some positive constant c’’, independent of n). 
Turning our attention back to the polyhedron Py» tet FIP.) = , 
U Fo(P.) denote the set of all faces of the polyhedron Pe The remainder 
of this proof will essentially consist of constructing a 1-1 mapping v 
from the elements of F(P,) into the elements of F(P?). We will then have 
that [F(P,)| < |F(P)]. 
Let Fi (P.JeF(P,) be any face of the polyhedron P,. If Fy l(P Lae 
(i.e. Fy(P.) is the empty face) then define P(F(P Jl=@. If, however, 
Fi(P.)#@ then let der Mtl be some point in the set F,(P,). Since 


UeF (PL) cP., we must have that cle; au; ; defines a legal capacity 
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function on the graph G, and that there exists some legal flow function f% 
on the network defined by G_ and cle; dau; js such that v(f“)eyro.,- Now 
2 
consider the point #’eR2" such that: 
Wi pRUy je for all 1si, jsn 
and 
7 ra . ' 
Witt fe; 5). for all lsi, jsn. 


>7 


Clearly W°eP” and thus there must exist some face of Pr, containing a’. 
Call this face Fx-(P7) and define P(Fy(P))=Fq-(P7). It should now be 
clear that W maps each face in F(P.) into a face in F(P7). Thus it 
_ remains to be shown only that W is a 1-1 mapping. Recall, however, that 
each face of a polyhedron is uniquely determined by the set of constraints 
its elements satisfy by equality. Further, notice that the points 4’ and 
U’ have the following relationship: 

1) wi je8 iff uf j=8, for all isi, jsn 

2) (AB iff 21348, for all 1sks2-2 

3) Wf j*8 for all {ij | e,,elm(t)} iff ypay=8. 
Thus for each distinct face F,(P.) of the polyhedron P_, WIFy(P,)) will be 
a distinct face of the polyhedron Re and so ¥ is in fact a 1-1 mapping. 
We can now conclude that |F(P,)| < |F (PA) | and thus by (3.6) we have: 
(3.7) |FIP,)| < glen?) 
Finally it should be clear that IF (Ph) I s |F(P,)|, for all s, and so by 


(3.7): 


272 
IF (PL ] s 20°), for alls. 
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3.5 Conclusions 

In this chapter we presented the linear decision tree model of 
computation [6,18,21,22], the notion of a polyhedron and a polyhedra! 
decision problem (21,22), and the class of polyhedral decision problems 
{B(P.) | n22} which most naturally arises when considering the maximum 
network flow problem. We then showed that the problem of establishing a 
non-linear lower bound on the linear decision tree complexity of max-fiouw 
can be reduced to the problem of establishing a non-linear lower bound on 
L(P,), for each n22. Next we demonstrated matching upper and lower bounds 
on the number of faces of dimension 1 composing each of the polyhedra P,, 
(n22). Finally, we established a 20(n?} upper bound on the number of 
faces, of any dimension, composing each of the polyhedra Pre Based on our 
results, we can now conclude that Theorem 3.1 can be of no use in 
establishing a non-linear Jower bound on max-flow, through the class of 
polyhedral decision problems {B(P,) | n22}. It remains an open question, 
however, whether or not the polyhedral technique can in genera! be useful 
for establishing non-linear lower bounds on max-flow. There may, for 
example, exist some more complex class of polyhedral decision problems that 


can be reduced to max-flou. 


{1} 


{2] 


[3] 


(4) 


(S] 


{6} 


[7] 


[8] 
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